稀醋酸/水溶液中单个醋酸与水分子的氢键结构

在醋酸/水体系的工业分离中,溶液中的氢键对分离效率有很大影响.本文采用两种第一性原理方法,即从头算分子动力学模拟(AIMD)和量子化学计算(QCC),对由单个醋酸和不同水分子所组成聚合体的氢键相互作用进行了研究,采用极化统一模型和自洽反应场模型计算得到了聚合体在水溶液中的热力学数据.从QCC计算的气相和水溶液中的聚合自由能表明六元环在两种状态下都为最优结构,热力学数据反映出的各种结构的相对稳定性与AIMD模拟的环分布符合得相当一致.研究表明,由于存在醋酸和水分子间的氢键作用,稀醋酸/水溶液中的醋酸分离要比在浓醋酸溶液中困难得多.     

二甲基精氨酸二甲胺水解酶-1的理论研究

运用分子对接和分子动力学方法研究二甲基精氨酸二甲胺水解酶-1(DDAH-1)与其抑制剂亚胺基烯丁基-L-鸟氨酸(L-VNIO)和亚胺基丙基-L-鸟氨酸(Me-L-NIO)的相互作用和结合模式,并根据实验得到的结论设计了亚胺基苯乙基-L-鸟氨酸(Ph-L-NIO)抑制剂.结果表明:L-VNIO比Me-L-NIO对DDAH-1的抑制效果更强,这个结果与实验测得L-VNIO和Me-L-NIO对DDAH-1的半抑制浓度IC50值大小一致.Phe75、Asp78、His172、Ser175和Asp268这五个氨基酸残基在三种抑制剂形成的复合物中起到非常重要的作用,从计算结果推断在这三个抑制剂中我们设计得到的Ph-L-NIO对DDAH-1的抑制效果最好.     

An N log N approximation based on the natural organization of biomolecules for speeding up the computation of long range interactions

Presented here is a method, the hierarchical charge partitioning (HCP) approximation, for speeding up computation of pairwise electrostatic interactions in biomolecular systems. The approximation is based on multiple levels of natural partitioning of biomolecular structures into a hierarchical set of its constituent structural components. The charge distribution in each component is systematically approximated by a small number of point charges, which, for the highest level component, are much fewer than the number of atoms in the component. For short distances from the point of interest, the HCP uses the full set of atomic charges available. For long‐distance interactions, the approximate charge distributions with smaller sets of charges are used instead. For a structure consisting of N charges, the computational cost of computing the pairwise interactions via the HCP scales as O(N log N), under assumptions about the structural organization of biomolecular structures generally consistent with reality. A proof‐of‐concept implementation of the HCP shows that for large structures it can lead to speed‐up factors of up to several orders of magnitude relative to the exact pairwise O(N²) all‐atom computation used as a reference. For structures with more than 2000–3000 atoms the relative accuracy of the HCP (relative root‐mean‐square force error per atom), approaches the accuracy of the particle mesh Ewald (PME) method with parameter settings typical for biomolecular simulations. When averaged over a set of 600 representative biomolecular structures, the relative accuracies of the two methods are roughly equal. The HCP is also significantly more accurate than the spherical cutoff method. The HCP has been implemented in the freely available nucleic acids builder (NAB) molecular dynamics (MD) package in Amber tools. A 10 ns simulation of a small protein indicates that the HCP based MD simulation is stable, and that it can be faster than the spherical cutoff method. A critical benefit of the HCP approximation is that it is algorithmically very simple, and unlike the PME, the HCP is straightforward to use with implicit solvent models. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Molecular dynamics simulation study of association in trifluoroethanol/water mixtures

Mixtures of Trifluoroethanol (TFE) and water with different proportions are studied using molecular dynamics simulations. The radial and spatial distribution functions, as well as the size distribution of TFE clusters are obtained from the trajectories. The variation of radial and spatial distribution functions with composition show that the addition of TFE enhances the water structure, but the hydrogen bonds between TFE molecules are broken as TFE is diluted with water. The TFE‐rich solutions have stronger TFE–water hydrogen bonds. The clustering of TFE molecules in low concentration region is attributed to the hydrophobic interactions between CF₃ groups. The distribution of cluster sizes in solution supports these conclusions. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Efficient nonbonded interactions for molecular dynamics on a graphics processing unit

We describe an algorithm for computing nonbonded interactions with cutoffs on a graphics processing unit. We have incorporated it into OpenMM, a library for performing molecular simulations on high‐performance computer architectures. We benchmark it on a variety of systems including boxes of water molecules, proteins in explicit solvent, a lipid bilayer, and proteins with implicit solvent. The results demonstrate that its performance scales linearly with the number of atoms over a wide range of system sizes, while being significantly faster than other published algorithms. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Polarizable empirical force field for sulfur‐containing compounds based on the classical Drude oscillator model

Condensed‐phase computational studies of molecules using molecular mechanics approaches require the use of force fields to describe the energetics of the systems as a function of structure. The advantage of polarizable force fields over nonpolarizable (or additive) models lies in their ability to vary their electronic distribution as a function of the environment. Toward development of a polarizable force field for biological molecules, parameters for a series of sulfur‐containing molecules are presented. Parameter optimization was performed to reproduce quantum mechanical and experimental data for gas phase properties including geometries, conformational energies, vibrational spectra, and dipole moments as well as for condensed phase properties such as heats of vaporization, molecular volumes, and free energies of hydration. Compounds in the training set include methanethiol, ethanethiol, propanethiol, ethyl methyl sulfide, and dimethyl disulfide. The molecular volumes and heats of vaporization are in good accordance with experimental values, with the polarizable model performing better than the CHARMM22 nonpolarizable force field. Improvements with the polarizable model were also obtained for molecular dipole moments and in the treatment of intermolecular interactions as a function of orientation, in part due to the presence of lone pairs and anisotropic atomic polarizability on the sulfur atoms. Significant advantage of the polarizable model was reflected in calculation of the dielectric constants, a property that CHARMM22 systematically underestimates. The ability of this polarizable model to accurately describe a range of gas and condensed phase properties paves the way for more accurate simulation studies of sulfur‐containing molecules including cysteine and methionine residues in proteins. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Mechanism of the decrease in catalytic activity of human cytochrome P450 2C9 polymorphic variants investigated by computational analysis

Cytochrome P450 (CYP) is deeply involved in the metabolism of chemicals including pharmaceuticals. Therefore, polymorphisms of this enzyme have been widely studied to avoid unfavorable side effects of drugs in chemotherapy. In this work, we performed computational analysis of the mechanism of the decrease in enzymatic activity for three typical polymorphisms in CYP 2C9 species: *2, *3, and *5. Based on the equilibrated structure obtained by molecular dynamics simulation, the volume of the binding pocket and the fluctuation of amino residues responsible for substrate holding were compared between the wild type and the three variants. Further docking simulation was carried out to evaluate the appropriateness of the binding pocket to accommodate substrate chemicals. Every polymorphic variant was suggested to be inferior to the wild type in enzymatic ability from the structural viewpoint. F‐G helices were obviously displaced outward in CYP2C9*2. Expansion of the binding pocket, especially the space near F′ helix, was remarkable in CYP2C9*3. Disappearance of the hydrogen bond between K helix and β4 loop was observed in CYP2C9*5. The reduction of catalytic activity of those variants can be explained from the deformation of the binding pocket and the consequent change in binding mode of substrate chemicals. The computational approach is effective for predicting the enzymatic activity of polymorphic variants of CYP. This prediction will be helpful for advanced drug design because calculations forecast unexpected change in drug efficacy for individuals. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Accurate free energy calculation along optimized paths

The path‐based methods of free energy calculation, such as thermodynamic integration and free energy perturbation, are simple in theory, but difficult in practice because in most cases smooth paths do not exist, especially for large molecules. In this article, we present a novel method to build the transition path of a peptide. We use harmonic potentials to restrain its nonhydrogen atom dihedrals in the initial state and set the equilibrium angles of the potentials as those in the final state. Through a series of steps of geometrical optimization, we can construct a smooth and short path from the initial state to the final state. This path can be used to calculate free energy difference. To validate this method, we apply it to a small 10‐ALA peptide and find that the calculated free energy changes in helix‐helix and helix‐hairpin transitions are both self‐convergent and cross‐convergent. We also calculate the free energy differences between different stable states of β‐hairpin trpzip2, and the results show that this method is more efficient than the conventional molecular dynamics method in accurate free energy calculation. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010